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ABSTRACT 

We present a detailed study of the dwarf galaxy VV124 (UGC4879), recently recognized as a remarkably isolated member of the 
Local Group. We have obtained deep (r^ 26.5) wide-field (23' x 23') g,r photometry of individual stars with the LBC camera at the 
Large Binocular Telescope under sub-arcsec seeing conditions,. The Color-Magnitude Diagram suggests that the stellar content of 
the galaxy is dominated by an old, metal-poor population, with a significant metallicity spread. A very clean detection of the RGB 
tip allows us to derive an accurate distance of D = 1.3 ± 0.1 Mpc. Combining surface photometry with star counts, we are able to 
trace the surface brightness profile of VV124 out to ~ 5' - 1.9 kpc radius (where fjr - 30 mag/arcsec^), showing that it is much 
more extended than previously believed. Moreover, the surface density map reveals the presence of two symmetric flattened wings 
emanating from the central elongated spheroid and aligned with its major axis, resembling a stellar disk seen nearly edge-on. We 
also present H i observations obtained with the Westerbork Synthesis Radio Telescope (WSRT), the first ever of this object. A total 
amount of ^ 10* Mq of Hi gas is detected in VV124. Compared to the total luminosity, this gives a value of Mm/Lv = 0.1 1, which is 
particularly low for isolated Local Group dwarfs. The spatial distribution of the gas does not correlate with the observed stellar wings. 
The systemic velocity of the Hi in the region superposed to the stellar main body of the galaxy is Vh = -25 km s"'. The velocity 
field shows substructures typical of galaxies of this size but no sign of rotation. The H i spectra indicates the presence of a two-phase 
interstellar medium, again typical of many dwarf galaxies. 

Key words. Galaxies: dwarf — Galaxies: Local Group — Galaxies: structure — Galaxies: stellar content — Galaxies: ISM — 
Galaxies: individual: UGC4879 



1. Introduction 

Until just a few years ago, VV124=UGC487SQ was con- 
sidered a unassuming isolated dwarf galaxy, classified as 
spheroidal/irregular, at a distance of D ~ 10 Mpc ( James et al. 
|2004 1. Integrated multi-color optical photometry was obtained 
by Jansen et al. ( 2004a) and Taylor et al. (|2005 [), and J and 



images were obtained by Grauer, Rieke & Quillen ( 2003| l, within 



large surveys of nearby galaxies. From inspection of a low res- 
olution integrated optical spectrum, Jansen et al.l (|2004b| con- 



cluded that it was likely "a young post-starburst galaxy". From 
the Ha +N[II] equivalent width, Jansen et al.| ( |2004bj l estimated 



Send ojfprint requests to: M. Bellazzini 
' A09125-I-5303 in the nomenclature adopted by 



(2004a I. 



Jansen et al. 



a total star formation rate of 0.005 M0yr ' (assuming distance 
of 10 Mpc). 

The generally adopted distance of D ~ 10 Mpc was based 
entirely on the re dshift estimate rep orted by the CfA survey 
{cz - 600 km s"'; Huchra et al. 1983 ). A team of russian sci- 
entists ( Kopylov et al.,2008, K08 hereafter), triggered by the ap- 
parent partial resolution of VV124 into stars in Sloan Digital 
Sky Survey images (SDSS, [Adehnan-McCarthy et aT][2^8) , 
carefully searched public databases and the literature, eventually 
making the case for a much lower recession velocity and a much 
smaller distance for VV124. They followed up this smart intu- 
ition with deep V,I photometry and low-resolution spectroscopy 
and were actually able (a) to resolve the galaxy into individual 
stars down to ~ 2 mag below the Red Giant Branch (RGB) tip, 
thus obtaining a direct distance estimate of D - 1.1 Mpc from 
the tip itself, and (b) to obtain a new estimate of the heliocentric 
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velocity, much lower than the CfA value, Vh = -70 + 15 km s"' 
(throughout the paper Vh stands for heliocentric radial veloc- 
ity). This meant that K08 had found a new member of the Local 
Group (LG), since with the newly determined distances and ve- 
locity VV124 is found to lie near the tum-around radius of the 
LG, and, in fact, being its most isolated member The galaxy 
has a remarkable elliptical shape and ranks among the bright- 
est LG dwarf spheroidal/transition type galaxies (Mb = -11.6). 
The stellar budget of the galaxy seems dominated by old stars 
(RGB; age ^ 2 Gyr) with colors compatible with low metal- 
licity (Z^ 0.001). However, a sprinkle of bright blue stars, and 
the identification of an Hii region, led K08 to conclude that 
VV124 is a transition type between dwarf irregulars (dirr) and 
dwarf spheroidals (dSph), like Phoenix, Antlia or LGS3 (Mateo 
|1998 ). The results by K08 were further discussed in more detail 
liTljikhonov et al. ( gUTO| TIO hereafter). 

According to K08, the location and the peculiar velocity 
of VV124 indicate that it has never been a satellite of a major 
galaxy of the LG, hence it evolved in full isolation for a Hubble 
time. Therefore, VV124 may contain a fossil record of precious 
information on the initial conditions of dwarf galaxies. It may be 
considered as a possible progenitor of the gas-less amorphous 
dSphs found in the vicinity of the MiUcy Way or M31, whose 
evolution has been likely largely driven by the strong interaction 
with the large galaxy they are orbiting around (see Mateo "1998'; 



Lokas et al. 2010a, and references therein). In particular, Mayer 
et al.| (|2006, 2007) have developed a detailed model, within a 
modern cosmological context, in which dSphs are produced by 
the morphological transformation of dwarf disk galaxies by tidal 
stirring and ram-pressure stripping during their path through the 
halo of the main galaxy they are gravitationally bound to (see 
also Kormendy et al.| [2008 and references therein for a more 
general view on the nature and origin of dSphs). As we shall 
see, VV124 may possibly share some remarkable characteristics 
with the precursors of modern dSphs envisaged in this model. 

The general interest in isolated galaxies as objects of undis- 
turbed evolution is witnessed, for example, by the large Hubble 
Space Telescope (HST) programme LCID ( |Gallart|2008[ ), aimed 
at the determination of the star formation history in the center of 
six isolated LG dwarfs of various morphological types. Here we 
are more interested in the structure and dynamics of a galaxy that 
should be untouched by the interactions with other large galax- 
ies since the beginning of time. In particular, the image presented 
in Fig. 1 of K08 suggests that the galaxy may be more extended 
than what could be enclosed into the 6' x6' field studied by those 
authors. For these reasons, we acquired much deeper observa- 
tions on a much wider field with the 2 x 8.4 m Large Binocular 
Telescope (LBT, Mt. Graham - AZ). A beautiful color image de- 
rived from these data is presented in Fig.[T] giving also an idea of 
the number and variety of background galaxies that can be found 
in such deep LBT images. In this paper we describe and discuss 
the results of these observations, as well as those from deep H i 
data obtained with the Westerbork Synthesis Radio Telescope 
(WSRT) and from low resolution optical spectroscopy obtained 
with the Telescopio Nazionale Galileo (TNG). 

The plan of the paper is the following: in Sect. |2]we present 
the LBT observations, we describe the reduction of these data 
and the artificial stars experiments. The process of surface pho- 
tometry of the innermost regions of the main body of the galaxy 
is described, and the adopted system of local coordinates is also 
introduced. In Sect.|3]we discuss the derived color-magnitude di- 
agrams (CMD), we provide a revised estimate of the distance to 
VV124 and we analyze the stellar content of the galaxy. Sect.|4]is 
devoted to the analysis of the surface brightness profile and the 



surface density distribution, while in Sect. |5] the results of the 
Hi observations are discussed in detail; the derived Hi velocity 
field is compared from the velocities obtained from low resolu- 
tion optical spectroscopy (Sect |5.l| l. Finally, the overall results 
are summarized and discussed in a broader context in Sect.|6] 

A few days before this manuscript was ready for submis- 
sion, a preprint was posted on the astro-ph archive (Jaco bs et al.| 
2010 hereafter JIO), presenting deep HST / Advanced Camera 
for Surveys (ACS) photometry of VV124. This study turns out 
to be complementary to ours, as it focuses on the star formation 
history (SFH) in the innermost ^ 40" of VV124, a region es- 
sentially out of reach of our photometry because of the extreme 
crowding (see Sect. 2.3 1. We will briefly refer to the results by 
JIO in the following, when appropriate, but we do not discuss 
them in detail. In general, for the issues treated in both papers, 
the results of the two studies are in good agreement. 

2. LBT Observations and data reduction 

Deep g and r photometry was acquired on the night of December 
2, 2008, at LBT, using the blue channel of the Large Binocular 
Camera ( jGiallongo et al. 2008 1. LBC optics feed a mosaic of 
four 4608 px x 2048 px CCDs, with a pixel scale of 0.225 ai'c- 
sec px^l. Each CCD chip covers a field of 17.3' x 7.7'. Chips 
1, 2, and 3 are flanking one another, being adjacent along their 
long sides; Chip 4 is placed perpendicular to this array, with its 
long side adjacent to the short sides of the other chips (see Fig. 4 
of Gianongo et al.|2()08 ). During our observations the pointing 
was chosen to place VV124 at the center of Chip 2, with the long 
side nearly aligned with the major axis of the galaxy (see Fig.[T]l. 
We will consider Chip 1 as our reference control field, sampling 
the back/foreground population in the direction of VV124; the 
CMDs from Chip 3 and Chip 4 are undistinguishable from that 
of Chip 1 , hence do not provide any additional information rele- 
vant for our purposes. For this reason they will not be discussed 
anymore in the following. In the following, we will use the terms 
Chip 1(2) and field 1(2), abbreviated as fl and f2, interchange- 
ably. Several long (f^.^^ - 300 s) and short (f^ jp - 20 s) exposures 
were acquired during the night, but we selected for the analysis 
only images taken under excellent seeing conditions, in particu- 
lar: four f„p - 300 s g band images with seeing ranging from 

300 s r band images with seeing 



0.60" to 0.69" and two t, 



exp 



ranging from 0.60" to 0.62". Three g and three r t^xp - 20 ex- 
posures were also reduced to provide a bridge between the pho- 
tometry from our long exposures (reaching the saturation level 
at r ^ 19.0) and the secondary calibrators from the SDSS R6 
catalog ( [Adelman-McCarthy et al.|2008| l that can be as faint as 
g,r ^ 23.0 but have average photometric errors crg,crr < 0.03 
mag only for g,r < 20.0. 

Relative photometry was performed independently on each 
image using the P SF-fitting code DAOPHOTII/ALLSTAR 
(Stets on|[l987[|199 4]l. Sources with peak higher than 3cr above 
the background were identified in a stacked image obtained by 
registering and co-adding all the images considered for the anal- 
ysis. Then, each of these stars was re-identified and fitted on each 
image (when possible). Only sources found at least in two g and 
two r images were retained in the final catalogue. The average 
and the standard error of the mean of the independent measures 
obtained from the different images were adopted as the final val- 
ues of the instrumental magnitude and of the uncertainty on the 
relative photometry. To clean the catalog from spurious sources, 
we also adopted cuts in the image quality parameters CHI and 
SHARP, provided by DAOPHOTII. After accurate inspection 
of the distribution of measured sources in the planes mag vs. 
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CHI and mag vs. SHARP, we decided to retain only sources 
having CHI < 2.0 and -0.5 < SHARP < 1.5. This selection 
removed 2893 (1185) sources, leaving a final f2(fl) catalog of 
15902 (7954) retained sources. 



2.1. Astrometry, Photometric Calibration and Reddening 

The instrumental positions, in pixels, were transformed into 
J2000 celestial coordinates by means of an astrometric solu- 
tion (in the form of a first degree polynomial) obtained with 
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CataXcon[^from 50(34) stars in common between the f2(l) cat- 
alog and the SDSS DR6 catalogue; the r.m.s. scatter of the solu- 
tion was 0.3"(0.2") in both RA and Dec. 

The instrumental magnitudes (gj, r,) were transformed into 
the SDSS ugriz absolute photometric system {gsoss, fsoss', g 
and r, in the following, for simplicity) with the following equa- 
tions 



gsDss = gi + O.lOSfe - n) + 34.359 
rsDss =r,-0.062fc-r,) + 33.911 



(1) 
(2) 



which were obtained from 43 stars with g < 20.0 in common 
between our f2 catalogue and the SDSS DR6 catalog; the rm.s. 
is 0.05 mag for both fits, implying an error on the zero point of 
^ 0.01 mag. The above equations have been obtained from stars 
covering the color range -0.37 < g - r < 1.5^ The calibration 
of the photometry of the other fields (1, 3, and 4) was obtained 
by fine adjustments (< 0.05 mag) of the zero point based on the 
handful of g < 20.0 stars in common with the SDSS per field, 
keeping fixed the color coefficients of Eq. [T] and |2] During the 
following analysis, in many instances we will have to compare 
observables in g,r with their counterparts in V,I. To do this we 
adopt 



V = g- Q.519ig - r) - 0.01 (cr = 0.005) 
B^g + 0.3l3(g-r) + 0.221 (cr = 0.011) 
from Lupton (200^, 

/ = r - 0.573(g - r) - 0.350 (cr = 0. 10) 
valid for g - r < 1 .4, and 



(3) 
(4) 



(5) 



g-r^ -1.417 + 2.650(y - /) - 0.600(y - if (cr = 0.06), (6) 

valid for 0.8 < ^ - r < 1.4. Eq. [5] and [6] have been derived by us 
from secondary standard stars in NGC 2419 (Bellazzini "2007). 

We interpolated the Schlegel et al.| ( [T998) reddening maps to 
obtain an estimate of E(B - V) for each source included in our 
final fl and f2 catalogues. For both fields we found an average 
E(B - V) - 0.015 and a standard deviation cte^b^v) - 0.001. In 
agreement with previous studies we conclude that in the consid- 
ered fields the reddening is very low and extremely uniform. In 
the following we always adopt E(B - V) - 0.015 and the red- 
dening l aws Af, ^ 3MEiB - V) and = 2.1\E{B - V), derived 
by ,Girardi et al. ( 2004) for cool metal-poor giants. 



2.2. Photometry of unresolved galaxies 

To get some characterization of the background unresolved 
galaxies which are so abundant in these deep and wide high 
Galactic latitude fields, we reduced the best g and the best 
r images with Sextractor (Bertin & Arnouts 1996| l; the as- 
trometric and photometric solutions were transferred from the 
DAOPHOTII catalog to this new catalog. Sextractor provides to- 
tal magnitudes, central surface brightness (SB), ellipticity, posi- 
tion angle and other useful parameters for extended objects in 



^ CataXcorr is a code aimed at cross-correlating catalogues and 
finding astrometric solutions, developed by P. MontegrifFo at INAF - 
Osservatorio Astronomico di Bologna, and successfully used by our 
group for the past 10 years. 

^ This range includes essentially all VV124 stars observed here (see 

* www. sdss .org/dr4/algorithins/sdssUBVRITransform. html 
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Fig. 2. Completeness fraction as a function of r magnitude for 
different color ranges (upper and lower panels) and for different 
(elliptical) radial ranges (different symbols). 



astronomical images. Moreover, it provides a "stellarity" index, 
based on neural network analysis, whose value runs from 5 = 1 .0 
for perfectly point-like sources (stars) to s = 0.0 for obviously 
extended sources (galaxies). By inspection of the g vs. s plot we 
found that the discrimination made by Sextractor between point- 
like and extended sources is safe down to ^ ~ 23.5 in our images, 
becoming increasingly blurred at fainter magnitudes. Hereafter 
we will refer to the catalog obtained with Sextractor and con- 
taining only sources with s < 0.5 as to the GAL sample, since it 
should be dominated by background galaxies at any magnitude 
and it should include almost exclusively galaxies for g < 23.5. 
The GAL sample will be used in Sect. [3] to interpret the back- 
ground contamination in the CMD of VV124. 

2.3. Artificial stars experiments 

The completeness of the stellar catalogs has been estimated 
by means of extensive artificial stars experiments. A total of 
100,000 artificial stars have been added to the images (follow- 
ing the recipe described in [Bellazzini et al."|2002a[ l and the en- 
tire data reduction process has been repeated as in the real case, 
also adopting the same selection criteria described above. The 
PSF adopted as the best-fit model for photometry was also as- 
sumed as the model for the artificial stars. Artificial stars were 
distributed uniformly in position, over the entire extent of f2, and 
in color, over the range -0.6 < g - r < 1.8 (see Fig. 5] below). 
They were distributed in magnitude according to a luminosity 
function similar to the observed one but monotonically increas- 
ing also beyond the limit of the photometry, down to r ^ 27.5 
(see Bellazzini et al. 2002a, for details and discussion). 

In Fig.|2]and Fig. [3] we show the completeness fraction (C/) 
as a function of r magnitude for different color ranges and for 



different (elliptical) radial ranges (see Sect. 2.4 for a definition 



of Tf), and the difference between input and output r magnitudes 
for different (elliptical) radial ranges, respectively. It is interest- 
ing to note that in the color range covered by the RGB of VV124 
(0.5 <g-r< 1.5, see Fig. |5]l, the completeness does not change 
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Fig. 3. Difference between input and output r magnitudes from 
artificial stars experiments in three different (concentric) ellip- 
tic annuli. The excess of stars with negative difference is due to 
blends: the recovered star is brighter because it is superposed to 
another source. 



with for > 1.5'; in this range C/ is fairly constant and larger 
than 80% for r < 26.0. Finally, we note also that for < 0.5' the 
completeness is lower than 50% virtually at any magnitude and 
the photometric errors due to crowding/blending becomes very 
significant. In agreement with these findings, in the observed cat- 
alogue there are only 129 stars in this radial range, nearly all 
brighter than r ^ 24.5. For these reasons, in the following we 
will limit our analysis to stars having > 0.5'. 

2.4. Surface photometry and coordinate system 
We performed surface photometry on the innerrnost 72" of the 



best g and r images using XVISTA0 (see |Lauer|l985 



Federici 



|et al.|[2007] for details on the code and on the adopted proce 
dure). The derived surface brightness profile will be discussed 
in Sect 



3.3 



below. For the purposes of the present section, 
we show in Fig. |4] the radial profiles of the ellipticity (e) and 
of the position angle obtained with XVISTA by fitting ellipses 
to the observed light distribution (PA, measured anti-clockwise 
from north towards east). Over the considered radial range, both 
quantities are remarkably constant with radius, independent of 
adopted passband. For this reason, we adopted the average val- 
ues (+ the associated standard deviation) as our final estimates 
of e and PA, in particular e = 0.44 + 0.04 and PA = 84.2° + 10.5°, 



in good agreement with the values reported by K08, |Jansen et al. 
( |2fe4b| ) and lTayloretal] ( 120051 ). 



To obtain the coordinates of the center of the smooth ellip- 
tical light distribution of the main body of the galaxy (i.e. not 
affected by the asymmetric distribution of the young stars), we 
obtained a smoothed intensity contour within the central ^ 1' 
from the r and g images; next we fitted these contours with el- 
lipses having e - 0.44 and PA - 84.2 deg, and we took the aver- 
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Fig. 4. Ellipticity and position angle as a function of distance 
from the galaxy center along the major axis, from XVISTA sur- 
face photometry. Red open squares are from the analysis of the 
r image, blue open triangles are from the g image. Each plotted 
point is the mean of the considered quantities over a radial range 
of ±10 px^ ±2.2" around the position of the point, the error- 
bars are the associated standard deviations (as in Federici et al.| 
^007 ). In the lower panel we plotted only points having uncer- 
tainties in PA lower than 30°. Continuous horizontal lines mark 
the average levels of the parameters considered in each panel; 
the corresponding values are also reported. 



age of the coordinates of the centers of the ellipses fitted to the g 
and r images as our final center coordinates (ao^^oj^ The g and r 
centers differ by 1 .3"; we took the average of the two as our final 
result. Our center is 11 .2" different (approximately in the E-N-E 
direction) from the coordinates provided by NED (and reported 
by K08 in their Tab. 1) and ~ 6" west from the center plotted 
by K08 in their Fig. 2. Tests on several images convinced us that 
our ffc^o provide the best approximation for the center of the 
overall light distribution, as shown in the lower panel of Fig.[T| 
above. Adopting these coordinates for the center, we convert to 
cartesian coordinates X,Y(m arcmin) projecting the equatorial 
coordinates of each star {a,6) on the plane of sky as in van de 
|Ven et al.| ( |2006| 



X/ = -(10800/;r)cos(5) sin(Q- - ao) 

Yt = (1 0800/71) [sin((J) cos(5()) - cos((5) sin(5o) cos(ck - ffo)] 

with Xi increasing toward west and Yi increasing toward north. 
These Xi, Yi coordinated were then rotated by 90° - PA into a 
X,Y system where the X axis is parallel to the major axis (and 
Y parallel to the minor axis) of the galaxy. Finally, we defined 



* The final estimates are the result of the iteration of the described 
process: first XVISTA has been run by adopting a center position esti- 
mated by eye, thus obtaining an estimate of e and PA, then the center 
was estimated by fitting ellipses and XVISTA was run again with the 
newly determined center position. The process converged after the first 
iteration, i.e. the e and PA profiles were essentially unchanged by the 
adoption of new coordinates of the center. 
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Table 1. Observed and derived parameters of VV124 



Parameter 




value 


notes 


ao 




09:16:03.18 


J2000 


^0 




+52:50:31.3 


J2000 


lo 




164.6619° 


Gal. long. 


bo 




42.8864° 


Gal. lat. 


SGL 




47.6121° 


Supergal. long. 


SGB 




-15.01 10° 


aupeigal. lat. 


E(B-V) 




A A 1 C 1 l\ r\f\ 1 

0.015 ± 0.001 




(m - M)o 




25.61 ± 0.13 




D 




1.3 + 0.1 Mpc 




1 arc sec 




0.3 pc 


conv. factor at D=1.3 Mpc 


6 




U.44 ± U.U4 




PA 




o4.Z ± lU.J 






T 1 1 
Zl.l 


+ 0.10 mag/arcsec^ 


a 


yUg(O) 


21.0 


± 0.10 mag/arcsec^ 




y"v(0) 


21.2 


±0.15 mag/arcsec^ 








1 T /I " 


Sersic scale radius 




17.5" 


Sersic scale radius 






41.3 


b 


Hr(Re < Re) 


22.5 


±0.15 mag/arcsec^ 




n, 




3y. / 


d 


gtot 




13. J ± U.IU 




^ tot 




n + n 10 

1 J .»J ^ »J. 1 v/ 




Blot 




13.8 ±0.15 




V,o, 




13.2 + 0.15 




My 




-12.4 + 0.20 




Ly 


8 


.2+[ « X 10« LyQ 


total V luminosity 


Vh 




-25 ± 5 km s-i 


e 






+ 17 km s"' 


galactocentric velocity^ 


^ stars 


l.( 


5 ± 0.3 X 10' Mq 


stellar mass^ 


Mm 




8.7 X 10^ Mq 


gaseous mass 



Notes, a Observed values, not corrected for extinction (note that the cor- 
rection is smaller than the uncertainty on SB in all the passbands con- 
sidered here). * R'^ is the elliptical radius containing half of the galaxy 
light.^ ^y{Rc < Re) is the SB within the ellipse of that radius; note that 
this is the quantity that must be used in Fundamental Plane plots, not 
the SB at R^, see Gadotti (2009). '' rj, is the radius of the circular area 
containing half of the galaxy light; the same value is consistently ob- 
tained from both g and r images. Mean systemic heliocentric radial 
veloc ity derived from Hi observations. ^ Adopting the solar motion 
from Schonrich, Binney & Dehnen||2010| ) and V^or = 220 km s '. * 
Adopting = 2.0. B an V magnitudes have been obtained from their 
g and r counterparts using equations |4]and|3] 



the elliptical distance from the center of the galaxy (or elliptical 
radius r^) as: 



X2 + 



that is equivalent to the major-axis radius. The X,Y coordinate 
system and will be always adopted in the following analysis. 

3. Stellar content and Distance 

3.1. The Color Magnitude Diagram 

The total CMDs of f2 and fl are compared in the upper panels 
of Fig. |5] Our new photometry reaches r - 26.5, a full 2 mag 
deeper than that of K08. The RGB of VV124 is the dominant 
feature in the f2 CMD, running from (g-r, r)^ (0.5, 26.5) to (g-r, 
r)^ (1.2, 22.5) and having no counterpart in fl. The comparison 



reveals that the prominent vertical plume lying to the blue of 
the RGB, around {g-r) - 0, is equally present in fl and f2, 
indicating that it is not due to a population that is characteristic 
ofVV124. 

In this context, it is worthwhile interpreting the features that 
can be ascribed to fore/background contaminating populations 
which can be identified in the right panel of Fig. [5] To obtain 
an idea of the foreground contamination by Galactic stars, we 
run a simulation with the TRILEGAL Galactic model ( |Girardi| 
eraL][2005| ). The synthetic TRILEGAL CMD for a FoV of the 



same area as fl and f2 is shown in the lower left panel of Fig.[5| 
the effects of photometric errors and completeness have been 
properly included. The model reproduces remarkably well the 
nearly vertical narrow plume observed in any field sig-r ~ 1.2. 
This is the well-known red plume of local M dwarf stars. The 
sparse diagonal band of stars going from (g-r,r)~ (1.1,26.0) to 
(g-r,r)~ (0.2, 19.0) is made of Main Sequence (MS) stars lying 
at different distances in the Galactic halo. Also this sequence is 
clearly identified in the CMDs of fl and f2, for r % 23.0, where 
it is not hidden by other populations. The feature may appear 
as remarkably narrow and well defined in the f2 CMD, but a 
simple inspection of the spatial distribution of the involved stars 
reveals they are clearly not correlated with VV124, confirming 
their Galactic origin. 

The above discussion unavoidably leads to the conclusion 
that for r ^ 23.5 and g-r % 1.2 the contaminating popula- 
tion is completely dominated by distant galaxies. This is a well 
known characteristic of deep photometric surveys (see, for ex- 
ample |Ibata et aL]|2007| their Fig. 6, in particular). This inter- 
pretation is strongly supported by the CMD of the GAL sample, 
shown in the lower right panel of Fig.|5j all the f 1 sources plotted 
here have been classified as galaxies by Sextractor (small dots) 
or by the SDSS pipeline (x symbols). They form two broad se- 
quences in color matching very well the ranges spanned by the 
well-known red and blue sequences of galaxies (see jBlanton &| 
Moustakas|2009[ for references and discussion). 

A minor population of young Main Sequence / Blue Loop 
stars is very likely hidden in the strong background contamina- 
tion, in this region of the CMD, since we see the bluest stars of 
this population emerging from the blob of contaminating galax- 
ies at g - r < -0.3 (see also JIO). K08 identified a small popula- 
tion of bright MS stars as young as a few tens of Myr, displaying 
an asymmetric distribution with respect to the smoothly elliptical 
shape of the overall galaxy (see Fig.[T]i. In particular, they found 
that the brightest blue stars and one H ii region are confined into 
a relatively narrow sheet in the southern half of the central 30" 
of the galaxy. Most of the g - r < -0.3 stars detected here in 
f2 are fainter than / ~ 24, thus they were out of reach of the 
K08 photometry. Figure |6] suggests that they are likely the faint 
counterpart of the bright blue stars discussed by K08: it is clear 
that they are associated with VV124 and the majority of them 
are confined in to the same strip located ~ 10 - 20" to the South 
of the galaxy center We used the synthetic CMD technique (see 



Cignoni & Tosi.20 10, for references) to reproduce the observed 



number of these young stars (limiting the analysis to r < 25.0), 
with the main aim of obtaining a rough estimate of the star for- 
mation rate (SFR) at recent epochs. The total mass of stars with 
age < 500 Myr is ~ 1.4 x 10^ Mq, assuming a Salpeter Initial 
Mass Function and a metallicity Z=0.001, i.e. ~ of the to- 
tal stellar mass of the galaxy (see Tab. [T]|. The resulting average 
SFR ~ 0.0003 Mq yr"' is in reasonable agreement with that ob- 
tained by James et aL] ([2004 ) from the integrated //„ flux, once 
rescaled to the correct distance of VV124 (1.3 Mpc instead of 
10.5 Mpc), i.e. SFR=0.00008 Mq yr \ given the considerable 
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Fig. 5. Upper panels: CMDs of the entire Chip 2 (left) and Chip 1 (right) fields. The dashed line marks the limiting magnitude 
reached in the V,I photometry by K08. The average r, g - r uncertainties are plotted as errorbars. The comparison between the two 
panels reveals the steep and wide RGB of VV124, tipping at r ^ 22.5, that dominates the f2 CMD while is completely lacking in 
the fl one. The red polygon enclosing the RGB of VV124 is used to select stars for the star counts described in Sect. [4] Lower left 
panel: predictions of the TRILEGAL Galactic model (Girardi et al. 2005) for a field of the same area as fl and f2 in the direction 
of VV124. The model includes only Galactic stars. The effects of photometric uncertainties have been properly added. Lower right 
panel: CMD of the sources classified as galaxies by Sextractor in fl (small dots). The blue x symbols are the sources brighter 
than r=20.5 classified as galaxies in the SDSS catalogue in the same area. The vertical lines approximately enclose the color range 
spanned by Red and Blue Sequence galaxies, according to Blanton & Moustakas (2009 1; see also Balogh et al. (2009[l. 



uncertainties involved. The agreement is good also with the SFH 
derived by J 10 in the same age bin. 

While the presence of these stars and their asymmetric distri- 
bution indicate recent activity in the galaxy (possibly connected 
to the asymmetric structure of the H i density and velocity field. 



see Sect|5]l, the associated star formation episode produced only 
a minor component of the overall stellar mix populating VV124. 
This young population is more abundant in the innermost region 
of the galaxy and is best characterized with the HST photometry 
by J 1 0, who concluded that it is quite rare at distances larger than 
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Fig. 6. Sources having r< 25.0 and g - r < -0.3 are overplotted 
(red open circles) on a LBC g image of VV124 zoomed on the 
center of the galaxy. The main body of the galaxy is enclosed 
within a square of side 1-3'. The orientation is the same as in 



40" from the center and that it formed in the last 500 Myr. The 
same is true for the very sparse population of candidate AGB 
stars brighter than the tip. For these reasons, in the following we 
will focus mainly on the RGB stars that are the visible part of 
the old stellar population which seems to dominate the galaxy. 



Fig. 7. Determination of the r magnitude of the RGB tip. Left 
panel: de-reddened CMD for the indicated elliptical annulus. 
The two parallel dashed lines enclose the RGB stars selected to 
derive the luminosity function. The pair of horizontal segments 
mark the level of the tip, the vertical segment marks the average 
color at the tip. Upper right panel: LF of the RGB. The sharp 
cut-off corresponding to the TRGB is marked by a vertical seg- 
ment. Lower right panel: Sobel filter response as a function of r 
magnitude. The peak corresponding to the TRGB is marked by 
a vertical segment. 



3.2. Distance from the tip of tlie RGB 



The tip of the red giant branch (TRGB; Lee, Freedman 



of age= 12 Gyr and 4 Gyr and spanning the full range of avail- 
able metallicity, we found that for TRGB colors (g-r)o < L4 the 
average magnitude of the tip is M™*-"^ = -3.03 with a standard 
deviation of just 0.04 mag. This is in excellent agreement with 

the predictions of the independent set of isochrones by Dotter 

& Madore|[T993) is widely recognized as one of the most |itaLl ( |2008| see BOS for some example in other passbands). The 
straightforward a nd general applicable distance indicators (see low scatter m M™"^* in the considered range of age and colors is 
IBell azzinT 2008 BOS hereafter, for references and discussion). 
In the present case, it is also the only standard candle accessi- 
ble with our data to estimate the distance of VV124. Figure [7] 
shows the very clean detection of the TRGB of VV124 we ob- 
tained from our data: the tip is unequivocally detected by the 
Sobel filter (see Lee, Freedman & Madore 1993 BOS) as an 
obvious cut-off in the luminosity function (LF) of the RGB at 
fTRCB _ 22.61 + 0.06, where the reported uncertainty is the half 



width at half maximum (HWHM) of the highest peak in the filter 
response (lower right panel of Fig.|7]). It is interesting to note that 
the uncertainty of the TRGB position is a factor of 2 lower than 
reported in TIO, owing to the much higher resolution/accuracy 
of our own photometry. The average color of the RGB at the tip 
is (g-r)o = L2S + 0.02. 

While an empirical calibration of the absolute magnitude of 
the TRGB in the SDSS passbands is still lacking, in BOS and 
|Bellazzini| ( |2009[ l it has been shown that all available theoret- 
ical models agree in predicting a very modest dependency of 
M™"^*, M™"^^, and Mf**^^ on metallicity and age, at least for 
relatively metal-poor models ([Fe/H] < -0.7; this is also the 
main reason which makes Mj'^^^ so appeahng as a standard can- 
dle, BOS). In particular, considering |Girardi et al.| ( [2004| l models 



confirmed also by BASTI ( Pietrinferni et al 2004) models, but in 
this case the average magnitude of the tip is ^ O.OS mag brighter 
than in other models (see |Bellazzini|2009] l. The latter difference 
is still within the uncertainty in the absolute zero point of the cal- 
ibration (see B ellazzini et al.|2 004a BOS and references therein). 
Adopting m;."'-'^ = -3.03 + 0.04, rj**^* = 22.61 ± 0.06 and in- 
cluding an uncertainty of 0. 1 mag for the overall zero point we 
obtain (m-M)o - 25.61 ±0.13; taking BASTI models we would 
have obtained (m - M)o = 25.72, instead. 

Another route to estimate the distance from the tip is to 
convert the observed magnitude and colors from ro, (g - r)o to 
Iq, (V - I)o using the relations reported in Sect. |2] and then to 
adopt the empirical calibration: 



j^TRGB 



O.OSO(y-/)^-0.194(y-/)o-3.939 +0.12 (7) 



derived by BOS from the original calibration as a function 
of [Fe/H] obtained in Bellazzini, Ferr aro & Pancino] ( |200I|l and 
revised in Bellazzini et al., (j2004a^. We obtain (V - /)^* = 
: 0.06 and, in turn, M/"^" = -4.05 ±0.12, that coupled 



0.9S1 



with/^'^'^^ = 21.53 ±0.10 gives (m-M)o = 25.5S ± 0.1 6, fully 
supporting the conclusions reached above, based on the |Girardi| 
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let al.| ( |2004"l l and [Potter et al.|P008] l models. As our final best 
estimate we adopt (m - M)q = 25.61 + 0.13, corresponding to 
D - 1.3 ± 0.1 Mpc; this value will be always adopted in the 
following. At this distance 1 arcmin corresponds to 378 pc. 

The main reason for the marginal difference with respect to 
the estimate by TIO (D - 1.1 ±0.1 Mpc) resides in the brighter 
level of the TRGB found by TIO: this may be due (a) to a larger 
impact of blending in their lower resolution photometry, and/or 
(b) to a systematic error in the absolute photometric calibration 
which should not aff'ect our photometry as we have the standards 
in field. Further support for the larger distance estimate obtained 
here is provided by the excellent agreement with the estimate by 

jidB 

3.3. Metallicity 

In Fig. |8] we compare the CMD of VV124 in different ellipti- 
cal annuli with a grid of RGB ridge lines of Galactic globu- 
lar clusters from the set by Clem et al. (2008 1, converted from 
the u'g/rfi/z' system to the ugriz one according to Tucker et al. 
( |2006| l. From blue to red, the RGB templates are for: M92, at 
[Fe/H]=-2.16 (all the metallicities are in the scale by Carretta 
|& Gratton|1997l CG97 hereafter), M3 ([Fe/H]=-1.34) and M71 
([Fe/H]=-0.70). The reddening and distance moduli adopted for 



the templates are taken from Ferraro et al. ( [1999 1. The compar 
ison of the observed RGB with GC templates is the standard 
way to obtain estimates of the metallicity of an old stellar pop- 
ulation from photometric data (see, e.g. Harris & Harris|[2000 
[Bellazzini, Gennari & Ferraro 2005| and references therein). It 
is known to provide valuable constraints on the overall metallic- 
ity distribution when the considered RGB is dominated by stars 
with age comparable with that of classical GCs (^8-10 Gy:[^. 
The SFH derived by J 10 from HST data indicates that the vast 
majority of VV124 stars are older than 10 Gyr even in the in- 
nermost region, thus supporting the idea that the analysis is well 
suited for the present case . 

The first general remarks that can be made based in inspect- 
ing Fig. |8] are : (a) the majority of RGB stars of VV124 are 
remarkably metal-poor, with [Fe/H] < -1.0, and (b) the ob- 
served color spread at any magnitude is significantly larger than 
expected from mere photometric errors, hence there should be 
a significant spread in metallicity (see below); both conclusions 
were drawn also by K08 and are confirmed by J 10. The distri- 
butions of RGB stars in the four panels seem to indicate that 
the innermost regions host a larger fraction of relatively metal- 
rich stars compared to the external ones, in agreement with the 
typical metallicity gradient observed in dSphs ( [Harbeck et al 



2001[ l. It may be hypothesized that the larger number of RGB 



redder than the M3 ridge line in the innermost annulus is due 
to a larger fraction of blending in this most crowded region (see 
Fig.|3]l. However, this cannot be the case, as blending spuriously 
increase the luminosity of stars, thus mimicking bluer colors and 
lower metallicities, not higher On the other hand, higher crowd- 
ing also implies larger photometric errors and, consequently, 
broader color distributions in the inner regions. Tests performed 
using synthetic stars from our set of artificial stars experiments 
suggests that it is very unlikely that the observed changes of the 

' JIO reports {m - M)o = 25.67 ± 0.04. Note that the reported error 
does not account for the systematic uncertainty in the zero point of the 
calibration, that is of order 0. 1 mag i Bellazzini et al. 2004a BOS). 

« See 

www . lorentcenter . nl/lc/2809/324/ f riday/Bellazzini . ppt 
for some recent tests on the reliability of the technique. 




Fig. 8. RGB fiducial ridge lines are super-imposed to the CMD 
of stars in four different elliptical annuli around the center of 
VV124. The ridge lines are labeled with the metallicity of the 
associated cluster (in the CG97 scale) in the upper left panel. 



color distribution of RGB stars with radius is due to mere obser- 
vational effects. 

We obtained a rough estimate of the metallicity of each star 
by linear interpolation on the grid of ridge line s of Fig. [8] as 



done, for example, in Harris & Harris (20001, Saviane et al 



( [2000| , and [Bellazzini et al.[ ( [2003| inb Umit any kind of con 
tamination as well as the impact of photometric uncertainties we 
considered only the stars having 22.7 < r < 24.0 and enclosed 
between the ridge lines of M92 and M7 1 . The average metallic- 
ity and the standard deviation in the outermost radial bin consid- 
ered in Fig.|8] i.e. the one less affected by crowding and photo- 
metric errors, are {[Fe/H) = -1.5 (which we adopt as the typi- 
cal metallicity of the galaxy, in the following) and c^^^^^j - 0.34 
dex. Once de-convolved from the dispersion due to observational 



effects (cr, 



synth 
[Fe/H] 



0.13 dex, estimated from artificial stars), this 



leads to an intrinsic one-sigma metallicity spread of ~ 0.3 dex, 
typica l of dSph galaxies ( [Mateo| 1998] [Tolstoy, Hill & Tosi et"aL 
[20091 1. 

3.3.1. An extended Horizontal Branch? 

JIO clearly detect an RC (and an old Red HB) approximately 
peaking at a magnitude corresponding to our ro ~ 26.5; they at- 
tribute part of this population to stars with age between 0.5 and 
1 Gyr. However, this is only a minor component of their SFH, 
hence the majority of the core-He-burning stars they observe 
must be associated with the dominant old population. We do not 
see any sign of an RC population, even in the less-crowded outer 
regions of VV124 (r^ > 1.5')- The LF of the stars in the color 
range expected to enclose such a feature is shown as a long- 
dashed thin line in the right panel of Fig. |9] and it shows no hint 
of the characteristic peak due to an RC population (compare, 
for example, with the case of Leo I, which hosts a significant 
intermediate-age population, Bellazzini et al.|2004b Smolcic et 
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Fig. 9. Left panel: de-reddened CMD of the stars lying in the 
elliptical annulus 1.5' < < 5.0', with the (RGB and HB) 
ridge lines of the GCs NGC2419 (dashed blue line) and Pal4 
(continuous red line), de-reddened and shifted to the distance 
of VV124, as derived here. The ridge lines has been derived by 
us from public SDSS data. The thin vertical lines enclose the 
typical color range of HB (-0.2 < (g - r)o < 0.4) and RGB 
(-0.4 < (g - r)o < 1.4) stars; the latter range is expected to 
enclose also the Red Clump of an intermediate-age population. 
Right panel: the LF of the stars in the HB color range in the con- 
sidered annulus (thick continuous line) is compared to the LF in 
the same color range from an equivalent area of fl (dotted line). 
The thin long-dashed line is the LF of the stars in the RGB color 
range. 



al.|2007| ). This is likely due to the sudden drop of completeness 
occurring in our sample in this magnitude range, although radial 
age/metallicty gradients may also play a role. Probably the RC is 
less conspicuous in the outer regions of VV124, if they are pop- 
ulated by older and more metal-poor stars than those inhabiting 
the center, on average. 

Along these lines, it is suggestive to note, in Fig. [9j that 
the HB ridge line s of the Galactic globular clusters NGC2419 
{[Fe/H] - -2.1, Harris] ( [T996 i; distance and reddening from 
Ripepi et al.| ( |2007| l) and Pal4 ([Fe/ H] = -1.4, [Koch & Cote 
( [201 0'); distance and reddening from Harris (1996)), nicely fall 
on top an over-density that has no counterpart in the CMD of 
the control field (fl). This is clearly shown in the right panel 
of Fig. [9] where the LF of VV124 stars in the color range 
-0.2 < ig - r)o < 0.4 is compared with stars in the same range 
from a field of equal area from fl. While the LF of the con- 
trol field rise gently up to r ~ 26 where it drops down to the 
limiting magnitude, the LF of the considered annulus around 
VV124 shows a sharp rise of star counts at r ^ 25.6, peaking 
at r ^ 25.50 



We tentatively interpret this feature as an extended old HB 
population associated with VV124. If this interpretation is cor- 
rect, the galaxy must host a sizable population of RR Lyrae: un- 
fortunately, the entire set of observations used in this paper was 
acquired in less than one hour, making it impossible to search 
for candidate variable stars. It must be recalled that all the re- 
sults presented in this sub-section concern stars within one mag- 
nitude from the limit of our photometry, i.e. a realm where the 
completeness is rapidly dropping and the photometric errors are 
quite large. Yet we feel that they are worth to be reported, in the 
hope they may trigger further investigations. 



4. Structure 



In Fig. 10 we show the azimuthally-averaged major-axis r-band 
surface brightness (SB) profile of VV124. The profile has been 
obtained by joining the surface photometry on concentric ellip- 
tical apertures obtained with XVISTA out to = 72" ^ 1.2' 
(see Sect. \2.4\ , with the surface density profile obtained from 
star counts in rj*"] which are equivalent to star counts on el- 
liptical annuli (see Fe derici et al.pOOT) for a discussion of the 
procedure and details). For star counts we used likely VV124 
members selecting candidate RGB stars enclosed in the polyg- 
onal box shown in the upper left panel of Fig. |5] Star counts 
are limited to the largest ellipse (with the parameters reported in 
Sect. |2.4| l which is completely enclosed within f2, i.e. only to 
RGB stars having < 6.0'. The two profiles overlap with two 
points in the region between ^ 1' and ^ 1.2': the overlap region 
was used to normalize the star counts profile, shifting it to the 
same scale of SB of the surface photometry profila^ The SB 
values (not corrected for extinction) as a function of are listed 
inTab.|2] 

Surface brightness profiles of VV124 were previously ob- 
tained (in several passbands) by Jansen et al. ( |2004a[ l and by 
Taylor et al. ( 2005| ), from surface photometry limited to < 
1.5'. Once transformed to the proper photometric band (using 
Eq. [3] and |4|l, both profiles are in good agreement with ours 
in the overlapping region. The novelty of Fig. 10 is that using 
star counts in the outer regions we were able, for the first time, 
to trace the SB profile out to ^ 5' - 1.9 kpc, reaching fir - 
30.0 mag/arcsec", thus demonstrating that the galaxy is much 
more extended than previously believed. It is interesting to note 
the a Sersic (1968 ) model of scale-length Rs = 0.29' = 110 pc 
and n = 1 .3 (i.e., a nearly exponential profile, in agreement with 
TIO) provides a good fit to the observed profile over the whole 
considered range 0.0' < < 5.0'. 

The main structural parameters we obtained from the ob- 
served profile are listed in Tab.[T] To estimate the integrated mag- 
nitude, we computed the total flux within the maximum elliptical 
aperture we reach with surface photometry (r^ - 72"); next we 
used the best-fit Sersic model to find the fraction of light outside 
that radius. It turned out that a correction of -0.1 mag should 
be applied to the magnitude within = 72" to account for the 
contribution of these external regions. By transforming our inte- 
grated g and r magnitudes into an integrated B magnitude (with 
Eq.|4]i we obtain B,o, = 13.8 + 0.25 in good agreement with pre- 
vious estimates from [Taylor et al.| ( [2005| reported also by K08) 



' All the results presented in this section are unchanged if 
completeness-corrected LFs are considered instead of the uncorrected 
ones shown in Fig.|9] 



I.e., keeping fixed the values of 6 and PA derived in Sect. 2.4 
For < r the star counts profile is strongly affected by 



radial 

variations in the completeness factor. The surface photometry profile 
was transformed onto the standard scale of g, r magnitudes using Eq.[T] 
and|2] 
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Table 2. Observed Surface Brightness profiles of VV124. 
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Fig. 10. Surface brightness profile of VV124 in r band obtained 
by joining surface photometry (for < 72"; Sect. 2.4 1 and 
star counts (for > 72") profiles. The continuous line is the 
Sersic model that best-fits the surface photometry. The level of 
the background (bkg, dotted line) and the associated uncertainty 
(dashed lines) are also reported. The error bars from surface pho- 
tometry are typically smaller than the dimension of the points. 



and Jansen et aL] ( |2004a) i, B,o, 
is My 



13.7. The absolute V magnitude 
12.4, corresponding to Ly - 8.2 x 10^ Ly,0- 



4.1. Density maps 

The two dimensional density maps shown in Fig. 



11 



are even 

more interesting than the extended profile shown above. Both 
maps are obtained by estimating the density from star counts on 
a fixed regular grid of nodes spaced by 0.2' in both directions, 
with an adaptive algorithm that adjusts the spatial resolution ac- 
cording to the local density. To describe the algorithm, it is more 

In 



11 



convenient to refer to the map in the lower panel of Fig. 
this case, likely RGB stars of VV124 are selected by requiring 
they lie within the selection box plotted in Fig. 12 similarly to 
what was done for the SB profile. Next, for each node of the 
grid we find the distance D50 of the 50th star in order of dis- 
tance from the node itself. Following this, the local density (in 
stars/arcmin^) is computed by dividing 50 by the area of the cir- 
cle of radius D50 (see Silverman 1986; Be llazzini et al.|[2002b| 
for discussion and another application of the same concept). The 
average density of the background, computed using the same se- 
lection over the whole f 1 , is then subtracted to this number and 
the total 1-cr error of the background-subtracted density is com- 
puted, assuming Poisson statistics. In the maps of Fig. 
density is expressed in units of cr above the background. 

The map in the upper panel has been obtained in the same 
way, but adopting the matched filter (MF) scheme ( |Rockosi"et 
|al.|2002| instead of a selection box. At each grid node the den- 
sity is computed by summing the weights of the 50 stars within 
Dso; the weights are assigned according to the star position in the 
CMD within the weight pattern shown in greyscale in Fig. [12] 
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22.64 ± 0.04 


22.96 ± 0.05 


sp 


0.56 


22.90 ± 0.04 


23.18 ±0.04 


sp 


0.63 


23.05 ± 0.02 


23.48 ± 0.04 


sp 


0.71 


23.18 ±0.04 


23.94 ± 0.05 


sp 


0.78 


23.29 ± 0.05 


23.99 ± 0.04 


sp 


0.86 


23.45 ± 0.03 


24.06 ± 0.03 


sp 


0.93 


23.60 ± 0.03 


24.12 ±0.03 


sp 


1.01 


23.81 ±0.09 


24.17 ±0.03 


sp 


1.08 


24.15 ±0.17 


24.34 ± 0.02 


sp 


1.16 


24.24 ± 0.25 


24.54 ± 0.02 


sp 


1.23 


24.28 ± 0.27 


24.79 ± 0.02 


sp 


1.30 


24.62 ± 0.06 


24.92 ± 0.06 


sc 


1.50 


24.87 ± 0.06 


25.17 ±0.06 


sc 


1.70 


25.22 ± 0.07 


25.52 ±0.07 


sc 


1.90 


25.72 ±0.08 


26.02 ± 0.08 


sc 


2.10 


26.29 ±0.10 


26.59 ±0.10 


sc 


2.30 


26.70 ±0.12 


27.00 ±0.12 


sc 


2.50 


27.18 ±0.14 


27.48 ±0.14 


sc 


2.70 


27.23 ±0.14 


27.53 ±0.14 


sc 


2.90 


27.67 ±0.17 


27.97 ±0.17 


sc 


3.10 


28.00 ± 0.20 


28.30 ±0.20 


sc 


3.30 


28.37 ±0.25 


28.67 ±0.25 


sc 


3.70 


28.98 ±0.35 


29.28 ±0.35 


sc 


3.90 


29.70 ±0.59 


30.00 ± 0.59 


sc 


4.10 


29.87 ±0.65 


30.17 ±0.65 


sc 


4.30 


30.05 ± 0.73 


30.35 ± 0.73 


sc 


5.30 


29.97 ±0.61 


30.27 ±0.61 


sc 


Notes, sp = surface photometry; sc 


= star counts. 





The weights are obtained as a map of ratios of the number of 
stars from f2 to the number of stars from fl in a given position 
on the CMD, normalized to have the maximum equal to unity. It 
is clear from Fig. 12 that the MF approach is more conservative, 
as the bulk of the signal is provided by the most likely genuine 
VV124 RGB members, most of the weight being assigned to 
stars brighter than r - 24.5. This is the reason the associated 
density map is less smooth than the one derived from the se- 
lection box, which equally weights all the candidate RGB stars 
down to r = 26.0. In any case, we want to stress that we have 
obtained density maps using various adaptive and non-adaptive 
(fixed radius) density estimators and in all cases we obtained 
density distributions with the same features shown in Fig. \11\ 

The most striking features of the density maps are two sym- 
metric, relatively thin wings emanating from the edges of the 
elliptical (bulge-likep^ main body around X 3' and reach- 
ing X ^ 8', approximately aligned along the major axis of the 
galaxy. These structures have no obvious counterpart in the other 
dwarf galaxies of the Local Group. The overall density distribu- 
tion recalls a disk galaxy with a prominent bulge, seen nearly 
edge-on. The galaxy is quite gas-poor as a whole and no hint 
of Hi is observed in the wings (see Sect. |5]l. Hence, the wings 
would be interpreted as the dry remnant of an ancient disk that 



Or pseudo-bulge-like (see, e.g.. 
therein). 



Gadotti 



2009 



and references 
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Fig. 11. Upper panel: density map obtained with the matched 
filter (MF) technique, adopting the pattern of weights shown in 
greyscale in Fig. 12 Lower panel: density map obtained using 

The 



12 



RGB stars selected to lie within the box shown in Fig. 
levels of grey range between 3cr to 8cr above the background, 
in steps of Icr, from the lightest to the darkest tone of grey. The 
density depressions near the center of the main body are due to 
the low degree of completeness in the very crowded innermost 
< 0.5' region. 




Fig. 12. Selection box (continuous line) and weight pattern (in 



greyscale) adopted for the density maps shown in Fig. 1 1 Higher 
weights corresponds to darker grey areas, going from 0. 1 (light- 
est grey) to 1.0 (darkest grey), with step 0.1. Weight - 0.0 is 
assigned to stars outside the grey area, hence these stars do not 
contribute to the density estimate. 



Fig. 13. Background-subtracted surface density profiles of the 
Eastern (blue open circles) and Western (red filled circles) wings 
of the density distribution shown in Fig.[TT[ from counts of RGB 
stars selected with the same box used for the density maps (see 
Fig. [T2) i. Upper panel: profiles along the Y direction. Lower 
panel: profiles along the \X\ direction; the adoption of the ab- 
solute value of X allows to compare the Eastern to the Western 
wing in the same scale. The dashed line is the best linear fit to 
the distribution of points from both wings. In all the cases only 
stars within |y| < 1.5' and 4' < \X\ < 8' have been considered. 



evolved undisturbed until the total consumption of the original 
gas it was made of. In this context, it is interesting to recall 
that TIO interpret the structure of VV124 as a superposition of 
disks of different thickness, the thickest being populated by old 
RGB stars, the thinnest by young Main Sequence and Red Super 
Giants. Alternatively, the wings may be identified with the in- 
ner regions of two tidal tails (see Mufioz, Majewski, & Johnston 
|2008; Klim entowski et aL]|2009| and references therein), but 
this explanation seems unlikely, given the extreme isolation of 
VV124 from any other mass distribution (Both the MW and M31 
are ^1.3 Mpc away, and the nearest dwarf is Leo A, at a distance 
of ^ 700kpc; see JIO). 

To further check the significance of these very unusual struc- 
tures, as well as their actual association to VV124, we obtained 
the surface density profiles of the two wings along the X and Y 
directions, using stars selected to lie within the RGB box shown 
in Fig 



12 



In both cases, the background has been estimated 
from fl, adopting the same selections and procedures, and sub- 
tracted to the f2 profiles. Since the major axis of the galaxy is 
inclined by just 5.8° with respect to the E-W direction, in the fol- 
lowing we will refer to the wing at positive X as to the Western 
wing and to that at negative X as to the Eastern wing. In Fig. [13] 
the Y (upper panel) and \X\ (lower panel) profiles of the Eastern 
and Western wings are compared. There are several conclusions 
that can be drawn from these plots: 

- The over-density in the considered regions is clear and sig- 
nificant in both profiles and in both wings. 
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- The profiles of the E and W wings are indistinguishable, both 
in the Y and in the \X\ direction. Both Y profiles have a clear 
peak near Y - 0, i.e. at main the plane of symmetry of the 
whole galaxy. 

- We do not see any obvious truncation of the wings in the X 
direction, but it is hard to draw any conclusion on this, given 
the very low surface brightness of the considered structures. 
Extrapolating the linear fit to the observed profiles (dashed 
line), S = is reached at \X\ ~ 12'. Note that the fit is not 
intended to bear any physical meaning, it was attempted just 
to perform the above extrapolation to zero surface density. 

From the results described above it can be concluded that the 
wings are clearly associated to VV124, as it is very hard to con- 
ceive how any kind of unrelated structure in the background can 
present such a degree of symmetry in the Western and Eastern 
sides of the galaxy, as well as such a degree of correlation with 
the major axis of the galaxy. 

From this analysis it is not possible to draw firm conclusions 
on the actual nature of the wings. However, it should be noted 
that the strong peaks at F = shown by the density profiles 
in the Y direction seems more typical of a disk than of a tidal 
stream. The kinematics of stars in the whole galaxy and in the 
wings will certainly provide very useful insight into the origin of 
the wings, e.g., by looking for the signature of coherent rotation 
about the minor axis (see, for example, the case of the isolated 
Tucana and Cetus dSph |Fraternah et al.|2009HLewis et al.|2007| 
respectively). 

To obtain further insight into the possible tidal origin of the 
wings, we performed simple plausibility tests by comparing the 
observed size of VV124 to its expected tidal radius for two basic 
scenarios, described below. The tidal radius r, is defined as the 
cut-off radius in the density distribution of a system of mass m 
imposed by the tides from a mass M lying at distance D. In this 
case we will use the formula 



2 / m \l 



(8) 



which is appropriate for a logarithmic potential (see 
[Bellazzini et al.|2004c| for discussion and references). 

1. Let us consider the tidal force exerted on VV124 by the 
Local Group as a whole, as if the LG was a single giant 
galaxy whose center is located at the barycenter of the LG, 
adopted to lie at D = 1.2 Mpc from VV124 ( |van den Bergh| 
[2000). This is clearly too simple a model for the mass dis- 
tribution of the LG, but it may provide a rough upper limit 
to the tidal force that may be felt by VV124. Taking the 
stellar mass listed in Tab. [T] as the total mass of VV124 
(a conservative choice, in this context) we find that for the 
range spanned by the most recent estimates of the LG mass 
(2 - 6 X 10'^ Mq, see (van der Marel & Guhathakurta|2008 



and references therein) the expected tidal radius ranges be- 
tween 8 kpc and 13 kpc. This is significantly larger than the 
maximum size of the main body of VV124 (r^ ^ 2 kpc), 
and also larger than the size of the wings themselves, which 
reach ^ 3 kpc. Hence, it seems unlikely that the structure 
of VV124 could have been shaped by the tidal interaction 
with the LG as a whole or by one of its major components, 
as both the MW and M3 1 are more distant and less massive 
than the ideal super-galaxy considered here. 
We now test the case of a close interaction with a dwarf 
galaxy. In the hierarchical paradigm, the accretion of small 




S'''16=":jO' 18=^20' 16"'10' -(6"'00' IS^SO" ;5'"-;0' 
'Ighl A!icension (;3000) 

Fig. 14. Total H i contours on top of an optical image (with image 
cuts putting in evidence only the densest central part of the stellar 
body of the galaxy). Contour levels are 5, 10, 20 and 50 x lO'^ 
cm"^. 



groups of dwarf galaxies by a major halo should not be un- 
common (see |D'Onghia et al.|2009| and references therein), 
thus such a scenario is not unlikely, in principle. We con- 
sider the nearest dwarf to VV124 , i.e. Leo A, a di rr galaxy 
with mass M ^ 8 x 10^ Mq ( |Brown et al.||2007 l). Solving 
Eq.[8]for D, we found that for Leo A to impose a tidal radius 
as small as 3 kpc to VV124, the two dwarfs should have been 
as close as ~ 10 kpc in the past. It is likely that from such a 
close encounter the two galaxies should have emerged either 



as a bound pair or with high peculiar velocities ( Sales et al. 
2007 ), none of which is currently observed. 

To follow up on these simple considerations, we plan to 
search for isolated dwarfs showing signs of tidal interactions in 
recent simulations of the evolution of LG-size groups, performed 
in a realistic cosmological context. However, this is clearly be- 
yond the scope of the present paper 



5. Neutral Hydrogen 

In order to study the possible presence of neutral hydrogen in 
VV124, the galaxy was observed with the WSRT on April 13, 
2008. A bandwidth of 10 MHz (coiTesponding to ~ 2000 km 
s"') was used which was covered with 2048 channels, resulting 
in a velocity resolution of about 2 km s ' after Hanning smooth- 
ing the spectra. The integration time was 12 hr. Calibration and 
analysis of the data were done following standard WSRT recipes 
using MIRIAD software ( |Sault, Te uben, & Wright 1995). Several 
datacubes were made by varying weighting and tapering. Lower 
resolution cubes, which, in principle, are more sensitive to faint, 
extended emission, do not show more H i compared to the full 
resolution datacube. The spatial resolution of this datacube is 
13 X 16 arcseconds and the noise level is 1.2 mJy beam The 
faintest detected emission is at the level of 5 x 10'^ cm"^. 
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Neutral hydrogen was clearly detected in VV124. Figure [l4| 

shows the distribution of Hi in relation to the optical galaxy. 

The total H i image was derived by smoothing the full-resolution 

data to 30 arcsec resolution and clipping this smoothed datacube 

at 3 cr to make a mask which was applied to the full resolution 

datacube. The Hi is fairly well centered on VV124, with the 

peak of the emission about 1 arcmin W of the center To the 

SE the H I extends to larger radius and seems to form a tail-like 

structure. S" 

o 

The flux integral of the Hi emission is 2.2 ±0.1 Jy km g 
s"' which, for a distance of 1.3 Mpc, gives a total Hi mass of 

a 

o 



8.7 + 0.4 X 10^ Mq. Given the optical luminosity of VV124, this 
implies that, although this galaxy contains a significant amount 
of H I, it is relatively gas poor. In a study of Local Group dwarf 
galaxies, Grcevich & Putman (12009) find that Local Group 



dwarfs that are, Uke VV124, about 1 Mpc distant from the Milky 
Way and Andromeda typically have Mui/Ly - 1, while for 
VV124 Mui/Lv = 0.11. Only the Cetus and Tucana dwaifs ai-e 
more gas poor (with Mm/Lv below 0.01). 

In Fig. 15 we show the velocity field of the H i emission. This 
velocity field was derived using the intensity weighted mean of 
the masked cube that was also used for deriving the total Hi 
emission. The kinematics of the Hi does not show any signa- 
ture of overall rotation over the optical body of VV124. Such 
kinematics is seen in most dwarfs with a dynamical mass below 
10*^ Mo (Be gum et al. 2008 ). The systemic velocity of the Hi 
coinciding with the optical body is Vh - -25 ± 4 km s"', quite 
different from Vh - -70 km s ' reported by K08; a detailed 
comparison between the H i velocity field and the results from 

The tail-like exten- 



5.1 



optical spectroscopy is reported in Sect, 
sion to the SE seems to have a velocity gradient towards more 
negative velocities. 

The origin of this Hi tail is not clear. A possibility would 
be that the gas is displaced from the main body by the com- 
bined action of stellar winds. The small gradient in velocity 
(5 - 10 km s"') and the small mass of gas in the tail (^ 10^ Mq) 
would require a kinetic energy of ~ 10^' erg, at most. Such an 
energy can be provided by about 10 OB stars within a few mil- 
lion years with an efliciency of transferring kinetic energy to the 
ISM of 10%. However, if the H i tail is caused by a galactic wind, 
the latter should be highly asymmetric given that the tail is seen 
only on one side and also there are no signs of double peaks in 
the line profiles that could point to an expanding shell. A sec- 
ond possibility is that of ram pressure stripping. This would jus- 
tify the head-tail morphology of the HI in VV124 and also the 
smooth and small gradient in the tail. The problem with this in- 
terpretation is the emptiness of the field around VV124, which 
forces us to assume that the ambient medium responsible for the 
stripping should be the intra-group medium associated with the 
Local Group. We can estimate the required density of the ambi- 
ent medium in this way ( [Grcevich & Putman|2009| l: 



"IGM 



™"™cm-^ 



(9) 



where cr jjj is the measured H i velocity dispersion, assuming 
that it is somehow representative of the potential (see below), 
hhi is the gas volume density, Vlg is the velocity of the galaxy 
with respect to the LG medium, Vlg - 72.6 km s"' (computed 
with Eq. 1 1 by Tully et aL]|2008[ l. Using the numbers relevant 
for VV124, we estimate that the density of the ambient medium 
should be higm ~ 5 x 10"'* cm"-*, which is probably too large at 
the distance of VV124 from the centre of the Local Group (e.g. 
[Anderson & Bregman||2010] l. For example, this value is larger 
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Fig. 15. Velocity field of the Hi emission of VV124 (colourscale, 
heliocentric velocity in units of km s ') with overlaid the con- 
tours of the total Hi (same contour levels as in Fig. 14 1. The 
superposed blue ellipses are centered on the optical center, has 
e = 0.44 and have the semi-major axis equal to the optical 
effective radius (inner ellipse) and a=2', corresponding to the 
- 26.0 mag/arcsec^ isophote. 



than the density required to produce the Magellanic Stream (at 
mere 50 kpc from the MW) in the model by (Mastropietro et 
[aLl|2005| l. Hence, while both scenarios outlined above do not 
provide a fully satisfactory explanation for the H i tail and gra- 
dient, the stellar wind hypothesis appears more likely, as a low- 
intensity recent star-formation episode is actually observed to 
occur near the center of VV124 (see Sect. [3] and JIO). 

It is worthwhile to note that the distribution of H i does not 
appear to con^elate with the wing-like substructures discussed 
in Sect. |4.1| above. Moreover, Fig. 



15 shows that virtually all 
the detected Hi is enclosed within the //^ - 26.0 mag/arcsec^ 
optical isophote, having semi -major axis ^ 2'. The ji,- - 
30.0 mag/arcsec^ optical isophote has semi-major axis ^ 4', 
hence the Hi contours do not exceed the optical body of the 
galaxy. 

The Hi spectra do indicate the presence of a two-phase in- 
terstellar medium. In the regions where the column density is 
below 2 X 10^" cm"^, the spectra show a single component with 



a velocity dispersion of, on average, 11 km s" (see Fig. 16 1. 
However, in the regions where the column density is above this 
value, the spectra show, in addition, a narrower component with 
a dispersion of about 4 km s This two-phase character of the 



ISM is seen in many dwarf galaxies ( 


Young & Loll 1996 1997 


Begum et al.|2006 Ryan- Weber et al. 


2008|l. 



A rough estimate can be made of the dynamical mass of 
VV124 by assuming a spherical and isotropic Hi distribution. 
In this case, the dynamical mass is given by Mdyn = Sr-mcr^ jG. 
For a velocity dispersion of 11 km s ' and a size of the Hi of 
'"hi — 475 pc this gives Mdyn = 6.6 x 10^ Mq. This would mean 
that Mdyn/Mbaiyon - 4 and || ^ 8. Such a baryon content is very 
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V„„(kms-) 

Fig. 16. Two H I profiles showing the presence of a two-phase 
ISM in VV124. The profile taken over a 1 arcmin x 1 arcmin 
region of the low column density on the eastern side of VV124 
only shows a broad component, while the profile taken over the 
region where the column density is above 2 x 10^*^ cm"^ shows 
the presence of an additional narrow component 



1 arcmin 



Fig. 17. Slit position for our DoLoReS observations, superposed 
to a LBT g image. The main targets are labeled. 



typica l for small isolated galaxies such as VV124 ( Begum et al. 
[20051 1. 

5.1. Low Resolution Spectroscopy 

Triggered by the apparent mismatch between the very accurate 
systemic velocity obtained from H i data (Vh - -25 km s"') and 
the estimate from low-resolution (R - A/AA ^ 1000) optical 
spectroscopy by K08 (Vh - -70 km s"'), we applied for Director 
Discretionary Time (DDT) at the TNG, to obtain a new indepen- 
dent estimate of the radial velocity from optical spectra. We were 
awarded ^ 3 hr to observe VV124 with the low -resolution spec- 
trograph DoLoRe^ 

The camera is equipped with a 2048 x 2048 px^ E2V 4240 
thinned, back-illuminated CCD. The field of view is 8.6' x 8.6' 
with a 0.252 arcsec/pix scale. Observations have been acquired 
trough a long-slit 1" wide, using the VHRV holographic grism, 
which provides spectra in the range 4752A < A < 6698A, with 
spectral resolution R - 1527 at 5725 A. The slit was positioned 
as shown in Fig. [TT] The main targets of our observations were 
(a) the bright star bl2, already observed by K08, who classified 
it as spectral type F51a, and (b) the central source CI, which ap- 
pears extended in our LBT images, hence it may be a (nuclear?) 
star cluster; CI was not observed by K08. Both sources have 
a total magnitude V ~ 18, but CI is expected to give a lower 
S/N spectrum because of its extended nature. In addition, the slit 
crosses the whole galaxy nearly along its major-axis, thus spec- 
tra of the integrated light can also be obtained. 

The observations have been performed during the night of 
March 10, 2010. We acquired four texp = 1800 s spectra 
with the above-described set-up and slit position. In addition, 
ArH-KrH-NeH-Hg and He lamp spectra have been acquired imme- 
diately after each scientific exposure, to secure the most reliable 
wavelength calibration. Finally, proper bias and flat-field images 
were also obtained. The spectra of each source were corrected 
for bias and flat-field, extracted and wavelength-calibrated, and 
finally stacked together into a final spectrum using standard 



http: //www. tng . iac . es/instruments/lrs/ 
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Fig. 18. Spectra of the observed sources in logarithmic scale and 
arbitrary units. 



IRAF tasks (see, for instance |Galleti et al.||2007| and refer- 
ences therein). The co-added spectra of CI and bl2 are shown 
in Fig. 



18 



together with the integrated spectra of the region be- 
tween the two sources. 

Several (relatively) strong absorption line^are recognized 
and labeled in the spectrum of bl2, which is remarkably similar 
to that shown by TIO in the upper panel of their Fig. 6. On the 
other hand, deep absorption is the only clearly recognized 



A problem with the subtraction of a very strong sky line around 
5800A seriously affected the final shape of the Na absorption line at 
5896A in all the spectra shown in Fig. Tsj for this reason this line has 
been always excluded from the analysis. 
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Table 3. Radial Velocities from DoLoRes spectra 



Source 








[kms-'] 




bl2 


-44± 18 


cross correlation 


CI 


-86 ± 20 


cross correlation 


int. light 


-6 ±30 


Ha emission 



feature in the spectrum of CI, possibly suggesting a young age 
for this candidate cluster. Finally, the spectra of the integrated 
light, in addition to a weak Hp absorption show an obvious //„ 
emission, indicating the presence of some hot gas, as also found 
by K08 and TIO. 

The heliocentric radial velocities of bl2 and C 1 have been 
obtained by cross-correlation with a set of ten template^^ us- 
ing the IRAF f xcorr task, following the procedure described in 
detail in |Galleti et al.| (|2007 ). As our final value of Vh for each 
source we took the average from the ten semi-independent es- 
timates. The standard deviation was adopted as the error The 
estimates were repeated using the centroid of the strongest lines 
and the results were found to agree with those obtained from 
cross-correlation within the uncertainties. The velocity of the 
integrated light spectrum was obtained only from the centroid 
of the Ha emission line. The zero point of our velocity scale 
was checked against the centroid of sky emission line and found 
to be accurate within ^ +10 km s"' (note, however that the 
bluest accessible sky line was 5460. 74A Hg; hence the stabil- 
ity of the wavelength calibration solution could not be checked 
in the bluest portion of the spectra). 

Before presenting our new velocity estimates, it is worth to 
recall the results of K08 and TIO. They found Vh = -90 ± 
15 km s"' and Vh = -82 ± 15 km s"' for the two supergiants bll 
and bl2, respectively, and Vh = -70 +15 km s"' for the overall 
integrated light (absorption lines). From the [OIII]5007A emis- 
sion in the integrated spectra of different regions they found Vh = 
-71 + 10 km s"' for their region a-slit2, Vh = -47 ± 15 km s"' 
for region b-sYiil, and Vh - -54 +15 km s"' for region a-slitl. 
Finally, they obtained two independent spectra of the H ii region 
they identified to the South-East of CI, finding Vh - -36 km s ' 
from slitl and Vh = -55 km s"' from slit2. TIO conclude that, 
while all the above estimates are lower (sometimes significantly) 
than the mean systemic velocity derived from the Hi obser- 
vations (Vh - -25 ± 5 km s"'), they fall in the km s"' to 
-50 km s ' range spanned by the overall Hi distribution asso- 



ciated to the galaxy. However, the inspection of the H 
field shown in Fig 



15 



velocity 
the radial velocities 



reveals that, in fact, 
of the neutral Hydrogen associated with VV124 range between 
Vh - -20 km s"' and Vh - -40 km s"'. Moreover, the sources 
having Vh 70 km s"' from TIO optical photometry are pro- 
jected on the side of the H i distribution where the velocity of the 
gas is anywhere larger than Vh - -27 km s"'. Hence the corre- 
lation between the various optical velocities by TIO and the Hi 
velocity gradient is quite poor 

Our results, summarized in Tab. [3] do not add very much 
to the puzzling scenario presented by TIO. The newly derived 
velocity of CI matches very well the Vh < -70 km s"' values 
that TIO find for bll, bl2 and the integrated light. On the other 
hand the velocity we obtain from the integrated //„ emission is 



The templates are high S/N spectra of M31 globular clusters with 
accurately known radial velocities, observed with the same set up 



I Galleti et al. 2007 1. As the typical integrated spectral type of GCs is 
F I Harris|1996 l, they are appropriate for the present application. 



compatible within Icr with the H i mean velocity. bl2 is the most 
interesting case as (a) it is the source for which we got the high- 
est S/N spectra and the richest one in terms of usable spectral 
lines, and (b) the target is in common with TIO, thus providing 
a direct comparison between the two analyses. While the S/N of 
our spectra and the TIO one are very similar, the spectral reso- 
lution of our set-up is higher by a factor of ~ 1.5. Indeed, our 
estimate of Vh - -44 ± 18 km s"' is at ^ Icr from the Hi one, 
while the estimate by TIO is at 3.6cr from that value. 

We conclude that the newly obtained Vh measures may help 
to reconcile optical and radio estimates of the systemic radial 
velocity of VV124, while still leaving reasons of concern (see 
Sect|6] for dicussion). 

6. Summary and Discussion 

We have presented the results of deep wide-field photometry, 
optical low-resolution spectroscopy and Hi observations of the 
dwarf galaxy VV124=UGC4879, recently recognized as lying 
in the outer fringes of the LG. The main results of our analysis 
can be summarized as follows: 

- While a sparse population of young stars (age:$ 500 Myr) is 
observed in the inner region of the galaxy, the dwarf is domi- 
nated by an old population. There are indications of the pres- 
ence of age/metallicity gradients, with older and more metal- 
poor stars being more prevalent at larger distances from the 
center of the galaxy. 

- We used a very clean detection of the RGB tip to obtain a 
new and more accurate distance estimate with respect to K08 
and TIO, D - 1.3 + 0.1 Mpc, in excellent agreement with 
the very recent estimate obtained by JIO from HST data. In 
spite of the larger distance, VV124 can be considered as a 
(present or future) member of the Local Group of galaxies, 
according to the criterion adopted by Mateoj ( |1998| see also 
the discussion and references on the method), and illustrated 
in Fig. [19] here. 

- Independently of the actual membership to the LG, it is con- 
firmed that VV124 is the most isolated nearby galaxy, very 
likely never disturbed by even a weak interaction with other 
dwarf or giant galaxies during its whole lifetime (see K08 
and JIO, for further discussion). For example, the relatively 
high recession velocity (V^^^+lOO km s ~"') of the gas-poor 
isolated dwarf Tucana ( Fraternali et al.|[2009) may indicate 
that it has been ejected to large distances by a three-body 
encounter, as envisaged by 'Sales et al. ( 2007| l. Such a sce- 
nario cannot be invoked in the case of VV124, which has a 
much smaller galactocentr ic velocity (V^=+17 km s see 
also .D'Onghia et al.|200"9l|Kazantzidis et al.|2010, for other 
scenarios for the transformation of dwarf galaxies and dis- 
cussion). 

- The main visible body of VV124 is remarkably elliptical 
(e ^ 0.44), regular and smooth, if the asymmetrically dis- 
tributed minor young component is ignored. We were able 
to follow the SB profile of the galaxy out to a distance never 
reached before, ~ 1 .9 kpc, demonstrating that the galaxy 
is much more extended than previously believed. 

- In addition, we obtained surface density maps of RGB 
stars revealing the presence of two relatively thin symmet- 
ric wings, emanating from the western and eastern edges of 
the inner elliptical body and aligned along its major axis, ex- 
tending out to ^ 3 kpc from the center of the galaxy. These 
low SB features {jir ^ 30mag/arcsec^) have no counterpart 
in other known galaxies of the LG (but see below). 
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Fig. 19. Heliocentric velocities versus the cosine of the angle be- 
tween a galaxy and the apex of the Sun motion with respect to 
the center-of-mass of the LG (see |Mateo|| 1998 ' and references 
therein), for galaxies within 2 Mpc from the barycenter of the 
LG (Dlc < 2 Mpc) in the catalogue of |Tully et al.| (|2008|, 
plus VV124 (data from the present paper). The parameters of 
the solar motion are also taken from Tully et al., while the po- 
sition of the barycenter is from |van den Bergh] (j2000). Filled 
squares are galaxies with Die ^ 10 Mpc, open squares have 
1.0 Mpc < Dlc < 1-5 Mpc, X symbols have Duj > 1.5 Mpc. 
VV124 has Dkj ^1.3 Mpc. The continuous line is the locus 
of rest with respect to the center-of-mass of the LG, the dotted 
lines are at +lcr and +2cr from that locus (cr = 60 km s"', from 
Mateo| 19*98 and in good agreement with what found here from 
galaxies having Dlc < 1 -3 Mpc). 



- Hi emission was detected for the first time by our high- 
sensitivity WSRT observations. The Hi total mass (~ 1 x 
10'' Mq) is small with respect to the overall stellar luminos- 
ity (mass), making VV124 more gas-poor than typical iso- 
lated dwarfs. The density peak of the gas is slightly offsett 
with respect of the optical center (by ~ 400 pc) and the dis- 
tribution shows a tail in the SE direction, also corresponding 
with the region with the most negative velocity. This may 
be indicative of some outflow/inflow process possibly tied 
to the most recent star formation event (see J 10), since the 
youngest stars are also preferentially found to the SE of the 
galaxy center 

- It is important to recall that the observed Hi contours do 
not extend beyond the optical body of the galaxy as traced 
by RGB star counts. Also, the H i distribution do not show 
any obvious correlation with the surface density stellar wings 
described above. 

- The velocity field of the H i does not show any sign of overall 
rotation over the main optical body of the galaxy, while a 
velocity gradient is observed in the tail-like structure in the 
SE side. The velocity gradient may indicate either infall or 
outflow of H I. The projected size of the tail and the projected 
difference in velocity with the main body give a timescale of 
several times 10^ yr. This timescale is very uncertain (due to 
projection effects), but could be connected to some elevated 



activity in the galaxy some time ago (see JIO). It could be 
an outflow due to the most recent episode of star formation, 
or the remnant of some H i inflow which triggered the star 
formation. 

- The gas associated with VV124 lies in the velocity range - 
20 km s"' > Vh ^ -40 km s"', with a mean systemic velocity 
of Vh = -25 + 4 km s"' . The interstellar medium is found in 
two phases, a narrower component associated with the inner 
regions (at higher column density), and broader component 
found over the whole body of the H i distribution, having a 
typical dispersion of 11 km s"^ 

- The correlation between the velocity of various sources as- 
sociated to VV124, as derived from optical spectra by us, 
K08 and TIO, and the Hi velocity field is quite poor. While 
the observed differences can be accommodated within the 
uncertainties, it is a bit disturbing that the large majority of 
optical estimates lie at Vh < -40 km s"', beyond the lower 
side of the range of H i velocities, in particular if one consid- 
ers that they are now taken from two independent sources, 
i.e. TIO and this work. This means that one should consider 
the possibility that there is a real difference between the sys- 
temic velocities of the stars and of the gas, due to the hypoth- 
esized gas flow suggested by the asymmetric structure of the 
Hi. Coming Keck-DEIMOS observations of RGB stars in 
the galaxy will hopefully settle this issue. 

- Fig. 20 shows that the structural parameters of VV124 are 
consistent with the scaling laws obeyed by dSphs and dirrs 
galaxies (see Kormendy et al.|20"08 Tolstoy, Hill & Tosi et 
al. 2009] and references therein). It is interesting to note that 
it lies at the upper envelope of the My -j-iy relation, i.e. it has 
the brightest SB for its total luminosity, and in an interme- 
diate position between dIrrs and dSphs in My - rn relation. 
According to Penarrubia, Navarro & McConnachie| ( 2008 1, 
the evolution in a relatively strong tidal field would have led 
to a larger scale radius, a lower luminosity and a lower sur- 
face density, thus driving VV124 toward the loci of "gen- 
uine" dSphs, in these planes. This is in qualitative agreement 
with the possibility that VV124 may be a precursor of mod- 
ern dSphs that did not enter in the interaction-driven evolu- 
tionary path that produced the latter class of dwarf galaxies. 

Probably the most intriguing features of VV124 revealed 
by the present study are the low SB stellar wings; in particu- 
lar if they are interpreted as an ancient stellar disk seen edge- 
on and considering the extreme isolation in which the galaxy 
evolved. As mentioned in the introduction, the idea that present- 
day dwarf spheroidals originated from dlrr galaxies that were de- 
prived from their gas at early times, as a result of tides and/or ram 
pressure stripping within the potential of a much larger galaxy or 
of a group/cluster, is not new (see Mateo |19"98 Tolstoy, Hill & 
ITosi et al.'2009 and references therein). The model recently de- 
veloped by Mayer et al. ( 2006 2007) explicitly postulates dwarf 
disk galaxies as the progenitors of dSphs (see [Kazantzidis et 
aLl|2010 for a recent and thorough discussion). Indeed, obser- 



vational evidence of the existence of the expected intermediate 
stages of such transformation process is beginning to emerge, 
from study of various environments (see, for example, the disky 
dEs identified in the Virgo cluste r by|Lisker et al.|2007| . In the 
context of the LG, [Saviane et al. (20101 have recently reviewed 



the evidence for the presence of the relics of a disk in NGC 205, 
suggesting that the process of transformation of the disk into 
a spheroid is currently ongoing in this satellite of M31. Two 
independent theoretical studies have shown that some relevant 
observational properties of the disrupting Sagittarius dSph and 
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Fig. 20. VV124 (black filled circle) compared to other LG dwarf 
galaxies brighter than My = -7.0 in the My vs. logarithm 
of the half-light radius (lower panel) and My vs. central sur- 
face brightness (upper panel) planes. Different symbols are used 
for different morphological types, according to the classification 
criteria by |Mateo| ( |1998| l. Absolute magnitudes and half-light 
radii are taken from from [Wa lker et al. (2010) and from Kalirai 
2010| l; when lacking in these sources they have been 



et al. 



drawn from [Mateo' ('1 998[l. Su rface brightness estimates have 



been taken from Mateofi 1998 1, |McConn achie & Irwin 



[Martin et a"L| ( 2009 ), Zuck er et al.| (jlOOl) , [Irwin et al. 



(2006 



(2008 



and Whiting et al. (1999) , with this priority. All the values of 
surface brightness are corrected for extinction. 



of the associated tidal stream can be more easily explained if 
a disk galaxy is adopted as the progenitor of Sgr (Penarrubia 
2010b[ ). NGC 205 and Sgr may sim- 



et al. 2010 Lokas et al 



ply represent two different stages on a similar evolutionary path 
driving the transformation of a similar low luminosity disky pro- 
genitor into a dSph via the interaction with the main galaxy. In 
this framework, VV124 would be akin to such disky progenitors 
but never entered the transformation path. Rather, it evolved pas- 
sively in isolation, thus preserving its disk intact until the present 
day. Assuming we are looking at a rotationally supported stellar 
system, the disk of VV124 would be relatively thick, but this is 
normally seen in the faintest dirrs in the Local Group and nearby 
clusters (Sanchez-Janssen et al. 2010). The thickness of the disk 
does not reflect environmental effects, but is rather the result of 
pressure support from either internal feedback from star forma- 
tion and/or the cosmic ultraviolet ionizing background becoming 
increasingly more important for the energy balance towards in- 
creasingly lower masses. Kaufmann, Wheeler & Bullock (2007) 
have shown how, for galaxies having circular velocities < 30 
km/s hosted in halos with typical spin parameters (A < 0.05), 
an effective temperature (i.e. thermal + turbulent) of the ISM of 
a few times 10^ K, or equivalently an ISM velocity dispersion 
of ~ 10 - 12 km/s, is suflicient to produce a substantially thick 
disk (with major-to-minor axis ratio in the range 0.2 - 0.4) as 
an equilibrium configuration. This is because in such systems 



the gas velocity dispersion is already close to the virial tempera- 
ture (for circular velocities below 30 km/s the virial temperature 
is < 5 X lO'* K), which forces the gas to acquire a high pres- 
sure scale height at equilibrium. Stars in such a system form out 
of a pressure supported turbulent gas disk and, being collision- 
less, have no way to dissipate such motions later on. Again, in 
the Local Group, at the lowest luminosity end of dIrrs, there are 
a number of examples of dwarfs with Hi velocity dispersions 
around 10 km/s (SagDig, Leo A and GR8 being some of these), 
and whose low rotation velocities (< 15 km/s) implies a halo 
circular velocity well below 30 km/s, consistent with the picture 
just outlined. 

Another very important issue is to understand how the galaxy 
became gas-poor in absence of stripping mechanisms due to the 
interaction with other galaxies. Photo-evaporation of the gas af- 
ter re-ionization ( Bullock, Kravtsov & Weinberg 200Q is a pos 



sibility certainly worth further investigation (see, e.g. Susa & 
Umemura [2004 1. Supernovae feedback might also play a role; 
recent cosmological simulations of the formation of a gas rich 
dwarf that are finally able to produce a realistic exponential disk 
with no bulge (Governato, Brook, Mayer et al. 2010) show that 
outflows at high redshift (z > 1) can remove more than 2/3 of 
the baryons even in dwarfs with circular velocities exceeding 
Vcirc = 30 km/s. Finally, given the fragility of such a low mass 
galaxy, we cannot exclude that the interaction with intergalactic 
gas clouds in the Local Group could have caused stripping of at 
least a fraction of an already diffuse, loosely bound interstellar 
medium. 

In this context it is very interesting to recall that also the Dark 
Matter (DM) halo expected to embed VV124 should be virtually 
untouched since the epoch of its collapse: the kinematics of the 
stars in the main body and in the wings of the galaxy should 
probe the mass profile of this pristine halo out to a remarkably 
large radius (~ 3 kpc), possibly opening a crucial window on the 
initial conditions of DM halos of this low-mass scale. 
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